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We investigate the influence of dissipation on one- and two-qubit rotations in coupled semicon- 
ductor quantum dots, using a (pseudo) spin-boson model with adiabatically varying parameters. 
For weak dissipation, we solve a master equation, compare with direct perturbation theory, and 
derive an expression for the 'fidelity loss' during a simple operation that adiabatically moves an 
electron between two coupled dots. We discuss the possibility of visualizing coherent quantum os- 
cillations in electron 'pump' currents, combining quantum adiabaticity and Coulomb blockade. In 
two-qubit spin-swap operations where the role of intermediate charge states has been discussed 
recently, we apply our formalism to calculate the fidelity loss due to charge tunneling between two 
dots. 
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I. INTRODUCTION 

The control of quantum superpositions in coupled 
electronic systems has been suggested as a possible 
way to realize quantum logic gates in semiconductor 
structures. An example are coupled quantum dots^"^, 
where Coulomb interactions between electrons can be 
exploited^"^ to define very small eflFective Hilbert spaces 
such as that of two tunnel-splitted ground states, sepa- 
rated by a large energy gap from the remaining many- 
particle states^"-^^. 

If parameters of the Hamiltonian (like the tunnel- 
coupling between dots) are slowly changed as a func- 
tion of time, adiabatic control of the state vector^^, swap 
operations'^, and the controlled transfer from an initial 
to a final state'^ become possible. In addition, by cou- 
pling such a system to external electronic reservoirs, one 
can pump electrons through the system in a controlled, 
adiabatic manner which in principle can serve as a 'read 
out' of the state vector in form of an electric current. 
In the simplest case, this can be achieved through a si- 
multaneous variation of two parameters as a function of 
time. 

Both adiabatic control of rotations or swaps, and co- 
herent pumping of electrons evidentally are very sensi- 
tive to decoherence. The coupling to dissipative degrees 
of freedom such as phonons disturbs the coherent time- 
evolution and therefore leads to a loss of control over 
the desired superposition. Decoherence and dissipation 
in quantum XOR gate operations have been discussed re- 
cently by Thorwart and Hanggi'^ in a powerful numerical 
scheme. They found that properties like gate fidelities are 
very sensitive to the dissipative bath coupling constant, 
but only weakly depend on temperature. 

In this article, we quantitatively investigate the role of 
dissipation for one- and two-qubit adiabatic 'rotations' of 
an electron between two coupled quantum dots. For the 
electron charge one-qubit, we discuss a 'quantum electron 
pump' for dots which are coupled to external leads (elec- 



tron reservoirs). In our scheme, decoherence is mainly 
due to absorption of bosons in excitations of the instan- 
taneous ground state which leads to an exponential tem- 
perature dependence. For a sinusoidal pulse, we use an 
exact solution and perturbation theory in the phonon 
coupling to predict how quantum mechanical oscillations 
between the dots can be made visible in a 'read out' elec- 
tronic current, similar to the experiment by Nakamura 
et al. in a superconducting Cooper pair box'^ and the 
recent electron spin resonance scheme for dots by Engel 
and Loss'''. 

First, we compare results for both weak and strong 
coupling to the bosons and derive an analytic expres- 
sion for the influence of the bath on the fidelity of the 
operation. Second, we investigate two-qubit swap oper- 
ations in spin-based two-electron quantum dots. These 
spin qubits sensitively depend on charge decoherence due 
to intermediate states where charge has tunneled between 
the dots. Piezoelectric phonons coupled to the electron 
charge incoherently mix states in the singlet sector and 
lead to a loss of fidelity of the swap operation. 

Adiabatic transfer and pumping of charges through 
small metallic islands or semiconductor quantum 
dots'^"^' has already been demonstrated experimentally. 
Furthermore, in the strong Coulomb blockade regime of 
coupled quantum dots, the non-adiabatic coupling to AC 
fields in photo-assisted transport^'^^"^^ has been estab- 
lished in experiments^^^^^ . In the opposite regime of 
weak Coulomb correlations, experiments in open dots^* 
have demonstrated the feasability of an 'adiabatic quan- 
tum electron pump'. These systems can be described as 
non-interacting mesoscopic scatterers^^"^'^. 

A combination of strong Coulomb blockade and the 
adiabatic control of the wave function in a triple quan- 
tum dot has been suggested recently'^''. Furthermore, 
in superconducting Josephson junction qubits'^^, adia- 
batic quantum computation with Cooper pairs'^^ and adi- 
abatic controUed-NOT gates for quantum computation 
have been proposed by Averin^^. 
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The decoherence properties of adiabatic one- and two- 
qubit operations are closely related to the dissipative 
Landau-Zener problem (see, e.g., the review^^ for further 
references). In time-dependent operations, decoherence 
rates in general become time-dependent themselves and 
therewith usually involve the whole (or a large part) of 
the spectrum of frequencies in the effective spectral den- 
sity p(uj) of the bosonic bath. As we show below, adia- 
batic pulses can be chosen such that qubit rotations run 
at a constant energy gap A to the excited state whereby 
dissipation is due to p{iv = ^/K), i.e., one boson fre- 
quency only. This can be used to extract p{uj) from the 
pump current by using a series of pulses with different 
A. 

A suppression of p{ijj) at certain frequencies oj = ojq 
has been predicted in free-standing 'phonon cavities' due 
to symmetry and geometrical confinement^^. Here, we 
demonstrate that a 'Rabi' rotation pulse tuned to con- 
stant energy difference A(i) — Hluq effectively 'switches 
off' the decoherence in such systems, at least within 
second order in the coupling constant. Consequently, 
this defines a one-dimensional 'decoherence-free mani- 
fold' (curve) on the adiabatic groundstate energy surface 
of the system, which might be of interest for adiabatic 
quantum computation schemes suggested recently""^. 

The paper is organized as follows: in section II, we dis- 
cuss one-qubit rotations, introduce the time-dependent 
spin-boson model, and describe the adiabatic transfer of 
an electron between quantum dots with and without dis- 
sipation. In section III, we derive a perturbative, analyt- 
ical expression for the one-qubit fidelity, compare with 
the strong-coupling case, and suggest a scheme to ex- 
tract quantum oscillations as a 'read out' electron pump 
current. In section IV, we turn to spin-qubit swaps and 
charge decoherence for the gate discussed recently by 
Schliemann, Loss, and MacDonald. Finally, section V 
is a short conclusion. 



II. ONE-QUBIT ROTATIONS 

Adiabatic transfer in a two-level system consists in the 
rotation of a chosen initial state |in) into a final state 
|out) by an adiabatic variation of system parameters. 
The simplest example is a spin ^ in a slowly rotating 
magnetic field. 

Using the analogy with adiabatic steering in atomic 
three-level systems, the possibility of adiabatic pumping 
of electrons through triple quantum dots has been sug- 
gested recently'^"', using a suitable 'design' of instanta- 
neous energies through two time-dependent tunnel cou- 
plings. Here, we consider the transfer of an electron 
through a double quantum dot via quasi-stationary adi- 
abatic eigenstates. This 'charge qubit'*. Fig. (1) left, is 
defined by two electron states \L) and \R) with a time- 
dependent energy difference e{t) = ^^(t) — e_R(t) and cou- 
pled by a tunnel matrix element Tc(t), as described by 



T(t) 




FIG. 1: Left: Double dot with time-dependent energy level 
difference e{t) = £l(J) — eR{t) and tunnel matrix element 
Tc{t), connected to electron reservoirs. Right: Surface of 
the lower energy eigenvalue £_ of the two-level Hamiltonian 
Hg'\t), Eq. (1). To adiabatically transfer an electron from 
the left to the right dot, e and Tc are varied as a function 
of time as in Eq. (2), corresponding to the curve on the e_ 
surface. 



the time-dependent Hamiltonian 

H^'\t) = '-^<^z+Ut)a,, (1) 

with a, := \L){L\ - \R){R\ and cr^ := \L){R\ + \R){L\. 
Experimental control of (constant) s and has been 
demonstrated in double quantum dots^'"^^^. 

The initial state is an additional electron in the left dot 
with an energy sl well below the chemical potential of 
the left lead, the final state is an additional electron in the 
right dot which then (keeping Tc — 0) is lifted above the 
chemical potential of the right lead. Such a transfer cycle 
of the open system (coupled to electron reservoirs) in the 
Coulomb blockade regime requires a Hilbert space 
spanned by three states |0), \L), and ji?)^'^". Here, the 
two basis states \L) = \N + 1,M) and \R) = \N,M + 1) 
describe one additional electron in the left (right) dot 
above a ground state |0) = \N,M) ('empty state'). 

For the remainder of this section, we only consider the 
first part of the transfer cycle, i.e., the dynamics of the 
double dot isolated from the leads, and turn to the full 
transfer cycle including tunneling to and from the leads 
in section III C. 



A. Coherent Adiabatic Transfer 

An adiabatic transfer from left to right requires the 
simultaneous change of at least two parameters such as 

e{t) = eo + ^1 cosfii 
Tc{t) = -r,exp[-(t-io)'M- (2) 

This corresponds to a change of e{t) with a simultaneous 
switching of the tunnel coupling Tc{t) between the dots. 
The precise form of the pulse, Eq.(2), is not important 
and has been chosen for convenience here. 

The instantaneous, hybridized eigenstates of the iso- 
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FIG. 2: Inversion (cr^) for transfer from left ({(t^) = 1) to 
right {{(Jz) = —1) in a two-lcvcl system, Eq. (1). Inset: time- 
dependent tunnel matrix element Tc(t) and energy splitting 
£{t), Eq. (2). Energies (times) are in units of the ampli- 
tude Tc {h/Tc) in Eq. (2); the other parameters are to = 25, 
n = 7r/(2to), r = 10, eo = 0.1, ei = -1. The adiabatic 
approximation Eq. (4) is shown as dotted line. 



lated coupled quantum dot are 
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|±) = ^[±2Te|L) + (AT£)|i?)] (3) 

N± := V4|7;|2 + (At£)2, A:=v/^M^^7^. 

The eigenvalues e± = ±5 A of the coupled system rep- 
resent two energy surfaces over the Tr-e plane, the lower 
of which (ground state) is shown in Fig. la. The pulse 
Eq. (2) corresponds to a curve on the e_ surface. The 
corresponding change of the 'inversion' {(Jz)t is obtained 
from a numerical integration of the (coherent) equation 
of motion for the density matrix of the system. The re- 
sult is shown in Fig. 2, together with the form of the 
pulses Eq. (2). 

In accordance with the adiabatic theorem, the initial 
groundstate \L) of the system is rotated into the instan- 
taneous superposition |— ) of \L) and \R), Eq. (3), if the 
rotation is 'slow', i.e., VI,t^^ ,1^"^ <C A/fi,. In this case, 
the time enters as a parameter into the state |— ) which 
is used to calculate the approximate expectation value 



(fT.)ad = -e(i)/A(t), 



(4) 



which excellently reproduces the overall form of the nu- 
merically obtained {(Jz)t- The form Eq. (4) corresponds 
to Crisp's solution for the adiabatic following of an atom 
in a near resonance light pulse^^'^^. 

The exact solution exhibits the expected Rabi oscilla- 
tions with frequency A{t)/h around the adiabatic value, 
which arc strongest when the tunnel coupling is fully 
switched on. Due to Landau-Zener tunneling from the 
adiabatic ground state |— ) to the excited state |-|-), there 
is always a finite albeit small probability Pl for the elec- 
tron to remain in the left dot, i.e., the excited state after 



the rotation. Pl can be made exponentially small for 
large enough level splitting A and slow pulses, but de- 
pends on the exact pulse shape s{t),Tc{t). 



B. Electron-Boson Coupling 

The adiabatic rotation discussed above is an idealiza- 
tion and valid only for a two-level system completely 
isolated from its environment. In semiconductor quan- 
tum dots, the coupling to low-energy bosonic excita- 
tions of the surrounding electron system, to photons, and 
to phonons leads to deviations from the coherent time- 
evolution. These excitations are characterized by an ef- 
fective spectral density 



Q 



■ Wq), 



(5) 



of modes Q with frequency cjq and constants gq cou- 
pling to the electron charge in the dot. In this paper, 
we only consider the coupling of phonons to the charge 
density which has been found to be dominant coupling 
mechanism in double quantum dots^'^. Other coupling 
mechanisms like inelastic spin-flip scattering^ ^ lead to ad- 
ditional dephasing channels, so that our results can only 
be considered as a lower bound on dephasing. On the 
other hand, the time-dependent spin-boson Hamiltonian 
{h= 1 throughout) 



(6) 



Q Q 

employed here (and derived in Appendix A) is of suffi- 
cient general form for our results being transferable to 
similar dissipative one-qubit rotations in other two-level 
systems. 

The function p{ui) can be calculated explicitely if one 
assumes sharply peaked electron densities of negligeable 
width around the dot centers r^/ji, having a distance d 
and zero extension into the z (growth) direction. For 
bulk piezo-electric phonons, one obtains^° 
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— sm — 

U) \UJdJ 



(7) 



where 5 is a dimensionless coupling constant and loj^ — 
c/d, with c the velocity of the phonon mode. In general, 
a finite extension I of the electron densities in lateral 
or growth direction leads to a form factor that cuts off 
phonons with frequencies w>l/c ('phonon bottleneck'). 

A microscopic determination of p{lo) would require 
(apart from details of the microscopic clcctron-phonon 
interaction potential) exact knowledge of the many-body 
electron density in the dots. We argue that the as- 
sumption of relatively sharply localized positions between 
which the additional electron tunnels is justified by the 
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strong intra-dot electron-electron repulsion^^ . The origin 
of the oscillatory form of Eq. (7)^° lies in the double-slit 
like scattering of the phonons (waves with wave vector 
Q) at the electron charge that is delocalized between the 
two dots. 

In the spin-boson problem, one often assumes a 
parametrized form of the spectral density'^'*, 



Ps{ui) := guj^e "c0[u)) 



(8) 



where Wc is a high-frequency cut-off. Notice that the 

ohmic case s = 1 describes the microscopic form ppz{uj), 
Eq. (7), in the limit ujd/^ ^ 0. However, phonons 
interacting with the dot through othc;r c;oiipling mech- 
anisms like the deformation potential coupling, surface 
acoustic waves^^, or Rayleigh-Lamb waves in confined 
geometries'^'^ , lead to other forms of the spectral density 
p{lo). Throughout the discussion of the one-qubit dynam- 
ics, we will assume the simple ohmic form, Eq.(8) with 
s = 1, and return to the form Eq. (7) in the discussion 
of the two-qubit in section IV. 



C. Equations of Motion in Presence of Dissipation 

The coupling to the bosonic bath introduces decoher- 
ence and in general (exceptions are discussed below) leads 
to a loss of fidelity of the adiabatic rotation, i.e. even for 
very slow and adiabatic pulses s{t),Tc{t), the transfer 
from the left to the right state remains imperfect with 
the final expectation value (ctz) considerably deviating 
from —1. During the rotation, the pure adiabatic state 
decays into a mixture. 

We solve for the reduced density matrix p{t) of the two- 
level system coupled to external electron reservoirs^' 
Different techniques can be applied for weak boson cou- 
pling (perturbation theory) and strong boson coupling 
(polaron transformations or path integral in 'NIBA' 
approximation^^'^''). In general, to obtain the solution 
for imc-dependent spin-boson problems even numerically 
is a non-trivial task"**^. 

From the Liouville-von Neumann equation, the diago- 
nal elements of p{t) are easily obtained as 



d_ 
dt 

d_ 
dt 



PLL{t) = -iT,(t)[pLR{t)- pRL{t)] (9) 

+ IL [1 - PLL{t) - PRR{t)] - iLPLhit) 

pRR{t) = -iTcit)[pRLit)- PLRit)] (10) 

+ 7fl [1 - PLL{t) - pRR{t)] - jRPRR{t), 



where we have included the coupling to the electron reser- 
voirs (see Fig. 1) with the tunnel rates (j — L, R) 



7i 



(11) 



where fj is the Fermi distribution in lead j and Tj the 
rate for tunneling between dot j and lead j. 



The equations for the off-diagonal elements PRL{t) = 
P*lr{^) ~ (-^|p(OI-^) '^^^ t)e obtained only approxima- 
tively by performing a perturbation expansion in either 
9 or Te: 



1. Bom-Markov Approximation 

In the first approach, one starts from the basis of the 
hybridized states |±) and considers the term V := \Aaz 
in the Hamiltonian Eq.(6) as a weak perturbation. If e 
and Tc were constant, one could easily define an interac- 
tion picture with respect to Hq := H — V and proceed 
in the standard way, i.e. second order perturbation the- 
ory and tracing out of the bosonic degrees of freedoms. 
However, when s{t) and Tc(t) become a function of time, 
the time-evolution of the unperturbed system in general 
can not be written down analytically which makes the 
evaluation of the electron-boson terms very tedious. 

Here, we use an adiabatic approximation by regard- 
ing the time t in e(t) and Tc{t) as a parameter for the 
derivation of the incoherent (electron-boson) part of the 
master equation for pLnit)- We neglect memory effects 
of the bosonic system and derive pLR{t) up to second or- 
der in V (Born-Markov approximation). The bosonic 
environment enters solely via the correlation function 
of the operator A, Eq. (6), in the interaction picture, 
K{t) = (i(i)i(O)) which yields 



K{t)= / (Lj p{ij)[nB{oj)e"^* + {l+nB{oj)) e- 
Jo 



(12) 



where nB{uj) = [e^'^ — 1] ^ is the Bose distribution at 
temperature T (/? = l/fc^T). The result is 



j^PLRit) = (ie{t)-^^^]pMt) 



(13) 



+ iTe(i) {pRR{t) - pLL{t)) 

- l{t) PLRit) + 7+(t) PLlH) - ^-{t)pRR{t). 

Here, the coefficients 7 and 7± (suppressing the param- 
eter t in the following) are defined as 
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7+ := 



7- : 



1 

J dt {e'^ + AT^ cos At) Re{K{t)}, (14) 

-| / dt {s {1-cos At) - i Asin At) K{t), 
A Jo 

-| / dt {s {1-COS At) -iAsin At) K*{t). 
A Jo 



The coupling to the bosonic bath introduces a dephasing 
rate 7 of the off-diagonal element plr, 



7 = 27rgp(A)coth('^^), 



(15) 



and additional terms (7±) in the coupling to the diag- 
onals. In the ohmic case s = 1 of the spectral density 
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p{A) an additional contribution to 7 linear in the tem- 
perature T appears which, however, is unphysical in the 
sense that for the microscopic ppz {uj 0) cx oj'^ this term 
does not appear, cf. Eq.(7) and Appendix B. Note that 
in order to be consistent within the Born-Markov approx- 
imation, the scattering rates in Eq. (14) have to fulfill^^ 
7(±) < keT. 

One should also keep in mind that within our adiabatic 

approximation, these rates are time-dependent through 
the time-dependence of s(t), Tc{t), and A = A{t) = 
\/s'^{t) + AT^{t). Furthermore, the imaginary parts of 
7± can be understood as a rcnormalization of the tun- 
neling rate Tg, as can be seen from Eq. (13). At finite 
temperatures, they have to be computed numerically. 



2. Polaron Transformation 

One can derive an equation for the off-diagonal pLR{t) 
in a second, alternative approach by performing a uni- 
tary polaron transformation of the original Hamiltonian. 
This method basically is a perturbation theory in the 
tunnel-coupling T,, and suitable for strong coupling to 
the bosonic bath. The unitary transformation 

H^H:= e'He-^ S =\a,Y.— («Q ' «-q) (16) 



Q 



is applied to the Hamiltonian, leading to 

^=^cT,+T,(i)[Xa++X+a_] +/fB, (17) 

where a+ = \L){R\, a- = \R){L\, and a constant c- 
number energy term in H has been dropped. The trans- 
formation Eq. (16) has removed the linear term \(JzA 
and instead introduced the product 



X = nQ£>Q ^ , £)q(0) :=e-"Q- 



(18) 



of unitary displacement operators for the bosonic modes 
Q into the tunneling term. The master equation for p^n 
is now derived in this transformed picture in lowest order 
perturbation theory in Tg. The result is 



PLR{t) = - [ dt'e'I^' 
Jo 



ds s(s) 



IL +1R 



C{t - t')pLR{t') 



+ iT,{t') {C{t - t')pLL{t') -C*{t- t')pRR{t')} [19) 

where C{t) is an equilibrium correlation function with 

respect to the bosonic bath (inverse temperature /3 = 
l/ksT, spectral density piuj), Eq. (5)), 



C{t) := {X{t)X^)B -- 
p{uj) 



-*(f) 



(20) 



f 

Jo 



duj- 



(1 



cos UJ 



t) coth 



I3lo 



+ i sin uit 



D. Dissipative Adiabatic Transfer 

In order to discuss the effect of dissipation on the adi- 
abatic transfer, wc now numerically evaluate and com- 
pare the equations of motion for the density matrix p(t) 
in both the perturbative and the polaron transforma- 
tion approach. Here, wc consider a closed two-level 
system without coupling to external electron reservoirs 
{Tl = Ffl = 0) and comment on the 'open' system case 
('quantum pump') below. The boson spectral density 
p{u)) = Ps=i{uj), Eq. (8), always has ohmic form. 

First, the strong coupling equations Eq. (9) and Eq. 
(19) can be condensed into a closed equation for {crz)t as 



) (21) 



|(^.>t = -J^dt'^[l±{aMf±it,t' 

f±{t, t') := ±2Te(i)Tc(t')Re {e=^'i'*' dM')C{t - t')} • 

This single integro-differential equation''* can be solved 
by standard numerical techniques. Results from both 
the perturbative approach, Eq.(13) and Eq.(9), and the 
polaron transformation, Eq.(21) are shown in Fig. (3) for 
the weak coupling regime {g = 0.01), where the same 
form of pulses Eq.(2) as in the coherent case. Fig. (2), 
has been chosen. 

We first checked in a separate calculation that the in- 
fluence of the imaginary parts Im{7±}, Eq. (Bl), on the 
solution of the master equation Eq. (13) is negligible. 
The perturbative inclusion of dissipation at low temper- 
atures leads to the expected, small change of the time- 
evolution of {(Tz), with (az) saturating slightly above the 
value {(Tz) = — 1 of the coherent case. This means that 
in the dissipative case, the onc-qubit rotation, i.e. the 
transfer of the electron from left to right, is never com- 
plete even if the rotation is performed adiabatically. We 
will derive an analytic expression for this loss of fidelty 
in presence of dissipation in section III. 

In comparison to the perturbative approach, the po- 
laron transformation approach breaks down at weak cou- 
plings g and low temperatures where it predicts inver- 
sions below —1 towards the end of the rotation. This 
breakdown if well-known from the static case (Tc, s 
fixed). On the other hand, the results from the polaron 
approach agree fairly well with the perturbative results 
at higher temperatures, as can be clearly infered from 
Fig. (3) which again is consistent with the spin-boson 
dynamics for fixed Tc and s. 



III. FIDELITY OF ONE-QUBIT SWAP 

The previous discussion has shown that the adiabatic 
transfer remains incomplete in the presence of dissipa- 
tion, i.e. (az) can considerably deviate from its non- 
dissipative value —1 in the perfectly adiabatic case where 
no Landau-Zener transitions occur. Our aim in this sec- 
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kBT=0.5 (PT) 
kBT=0.5 (pert.) 
kBT=3.0 (PT) 
kBT=3.0 (pert.) 




FIG. 3: Two-level system as in Fig. (2), but with (weak) 
dissipation at finite temperatures T. Comparison between the 
strong coupling polaron transformation (PT) solution and the 
perturbative solution is made. 



tion is to analyze how the deviation 



(22) 



depends on temperature, boson spectral density p{lo)^ 
and the duration of the swap operation. Intuitively, it is 
clear that for too long swap duration, inelastic transitions 
to the excited level will have sufficient time to destroy the 
coherent transfer. On the other hand, if the swap opera- 
tion proceeds too fast, Landau-Zener transitions become 
stronger and again lead to deviations 5{(Jz) > 0. This 
means that for any given inelastic rate Fin there should 
be an optimal swap duration such that these two com- 
peting processes balance each other. 

One problem in quantifying this argument is that Fin 
naturally is time-dependent itself. Silvcstrini and Stodol- 
sky nevertheless have suggested**^ to extract decoherence 
rates in experiments by performing an adiabatic sweep 
and determining the maximum sweep duration for which 
the inversion is still successfuU. However, their model of a 
Bloch equation (with Fin as a constant parameter) is too 
simple in order to make detailed quantitative predictions 
here. 

In the following, we calculate the deviation 5{az)t for 
both weak and strong electron-boson coupling and derive 
analytic expressions that allow to quantify the arguments 
above. The main physical idea in this discussion is the 
introduction of curves (ellipses in the e-Tc plane) where 
the excitation energies 



A{t) = ^/e{tf+An{t)^ = A 



(23) 



remain constant as a function of time and consequently, 
dissipation is determined by a constant boson spectral 
density /9(A). 



A. Transfer Rate: Exact Solution for Rabi 
Rotation and Weak Dissipation 

In the weak coupling case, wc follow Grifoni and 
Hanggi^* and perform a unitary transformation of the 
original Hamiltonian (6) to a Hamiltonian H := UHU~^, 



^^'T""^" V2A^^+ A^^ 



A + Hb, 



(24) 



hy- 



where the matrix U contains the columns of the 
bridized eigenstates |— ) and |-|-), Eq. (3), dz = |— )(— | — 
I +)(+!, cFx = |— )(+| + |+)(— I, and the time-dependence 
of Tc{t) and £{t) is again considered as parametric. We 
introduce the amplitude a_| (t) for an inelastic transi- 
tion from the adiabatic ground state |— ) at time t = to 
the excited state at time f > 0, i. e. a transition in- 
duced by the coupling to the bosons. The corresponding 

probability {t) = (t)p = ■^S{az)i describes the 

deviation 5{az)t of the inversion due to the coupling to 
the bosons. The amplitude a-\ (t) is given by the non- 
vanishing matrix clement of the time evolution operator, 
expanded to lowest order. 



o+_(i) = i [ dt' 
Jo 



A(i') 



(25) 



where the interaction picture is with respect to —A/2az + 

Hb. With {+\ax{t')\-) = ex.p(ij*' dsA(s)), we find the 

corresponding probability for a transition from — to -|- 
due to the interaction with the bosons as^^ 



= j dwp(a;)|nB(w)/(w,t) 
+ [l + nB{uj)]f{-iJ,t)] 



f dt' 
Jo 



A{t')' 



-iJo ds[A{s)-u;] 



(26) 



'Elliptic' pidses {Tc{t), e{t)), Eq. (23), arc defined by 
curves on the £_-surface (Fig. 1) with constant energy 
difference A to the excited state e+. For the particular 
sinusoidal form 



Tc{t) 



■smflt, s{t) = -Acosnt (27) 



of Tc{t) and £{t), the time-dependent Hamiltonian 
ijW(^) = -(A/2)[cos(m)£7^ -I- sin(m)(7^], Eq.(l), be- 
comes exactly integrable. It corresponds to a spin ^ 
in a magnetic field that rotates within the .x-z-plane 
around the y-axis with frequency f2. The solution for 
the inversion {az)t is easily obtained by transforming the 
Schrodinger equation into a rotating frame (Rabi solu- 
tion). One obtains 



Rabi 



+ cos LORt 



cosnt (28) 



-|- — sinwijf sinOf, w/j := -\/fP"-i-A^. 
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prediction Eq.(31). The agreement between the numer- 
ical and the analytical prediction (dotted lines), is ex- 
tremely good for small coupling constants {g = 10^^), 
but deviations become strong at small Q where the simple 
second order perturbation theory becomes worse. The 
1/n dependence of the dissipative contribution to S{az) / 
is clearly visible at small ^, indicating that for too long 
pulse duration the electron swap remains incomplete due 
to incoherent dissipation. On the other hand, if the pulse 
duration is too short (larger Q), the oscillatory coherent 
contribution from (fr^)^''^' dominates. 

One should bear in mind, however, that Eq. (31) is an 
approximation that only holds in the limit of an infinitely 
slow adiabatic change, i.e. tf = n/fl oo. In fact, for 
any finite pulse duration tf < oo, even in the limit of 
zero temperature T = 0, Eq. (26) yields 



FIG. 4: Inversion change 5{az)f after time tf = n/Q for sinu- 
soidal pulses Eq.(27) as obtained from the master equation, 
Eq.(13). Dotted curves correspond to the analytical predic- 
tion, Eq. (31). 



Here, the Rabi frequency luu corresponds to maximal de- 
tuning (5 = f2 — since there is no static 'magnetic field' 
oc ay in ^/-direction. 

The use of the harmonic pulse Eq. (27) has the further 
advantage that the quantity f{oj,t), Eq. (26), can be 
evaluated analytically. The swap operation requires a 
pTilsc acting half a period from the initial time i = to 
the final time tf = 7r/f2. Using t — tf in f{(jj,t), we 
obtain 



P+-{tf) 



f 

Jo 



(L)p{ui)f{-ui,tf), 



(32) 



which shows that there is a small, but finite probability 
for inelastic processes during transitions from |— ) to |-|-) 
even at zero temperature, in agreement with^^. These 
processes are due to the spontaneous emission of bosons 
which occur during Landau- Zener transitions from |— ) to 
|-|-) with a finite probability as long as t/ is finite. 



In the adiabatic limit f2/A 0, we find an approxima- 
tion to the integrals Eq. (26). In that limit, one has 



whence the dissipation induced change of the inversion 
becomes S{a,)f'''' = 2P+_(t/) = 2f^p{A)nB{A). The 
sum of the coherent contribution 6{az)t from the Rabi 
solution at t = tf, Eq. (28), and the perturbative dissi- 
pative contribution 5{az)'f^^^ leads to 



B. Transfer Rate: Crossover for Strong Coupling 

For strong electron-boson coupling constants g, we 
used Eq. (21) to extract the final inversion (ctz)/ af- 
X := {A- u;)/fl. (29) ter the time t = tf for the pulse Eq.(2). An interesting 

observation can be made in the temperature behavior of 
{(7z) f. for couplings < 2, a temperature increase leads 
to an increase of {cTz) f, which is as in the weak coupling 
case. However, above g > 2, the temperature dependence 
changes in that larger temperatures T lead to smaller val- 
ues of {<Jz)f- In fact, for large coupling constants g, the 
system tends to remain localized in the left dot state \L) 
and no timneling to the right state \R) occurs. In this 
regime, higher temperatures destroy the localization and 
lead to smaller {(Jz)f- This behavior again is consistent 
with the borderline g = 2 (a = 1)'*^ in the dissipative 
two- level dynamics^'' for static parameters e and T^. 



4\/2 



(30) 



AV + fiLVcos(^) 



(31) 



+ 2 



P(A) 



ne^p{A/kBT) - 1 



n<^A, c = 2.4674. 



Results for the; inversion change S{(7z) y from the ideal 
value (az) / = — 1 as a function of the pulse frequency fl 
for a half-period sweep (duration tf = ir/fl) are shown 
in Fig. (4). We compare the results from the (weak- 
coupling) master equation, Eq.(13), and the analytical 



Adiabatic Quantum Pumping in Open Double 
Dots 



In the remainder of this section, we discuss how results 
for the closed double dot system, and in particular the 
coherent and incoherent inversion swap deviation S{az) f 
like the one for the Rabi rotation, Eq. (31), relate to 
electron transport in a double dot system coupled to ex- 
ternal leads. Here, we only consider the weak dissipation 
case. 
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FIG. 5: Inversion (uz) for strong electron-boson coupling after 
application of the pulse Tc{t) and e(t), Eq. (2), as in Figs. 
(2) and (3). Clearly visible is the crossover at g ~ 2 where 
the temperature dependence changes. 




FIG. 6: Adiabatic scheme for cycling electrons from left to 
right through a double quantum dot with energy level differ- 
ence e{t) = £L{t) — eR{t). Left and right energy levels are 
assumed to change symmetrically with respect to e = 0. 



The main idea is to apply time-dependent pulses such 
that the quantum mechanical time evolution of the two- 
level system is well separated from a merely 'classical' 
decharging and charging process. One complete cycle of 
such an operation is sketched in Fig. 6. The cycle starts 
with an additional electron in the left dot and an adia- 
batic rotation of the parameters {e{t), Tc{t)) such as, e.g., 
in Eq. (27) or in Eq. (2), cf. Fig. 6 a-b. This completely 
quantum-mechanical part of the cycle is performed in the 
'save haven' of the Coulomb- and the Pauli-blockade^^, 
i.e., with the left and right energy levels of the two dots 
well below the chemical potentials /i of the leads (which 
are assumed identical here for simplicity) . The cycle con- 
tinues with closed tunnel barrier Tc = and increasing 
eii{t) (Fig. 6 c); the two dots then are still in a super- 
position of the left and the right state. The subsequent 
lifting of the right level above the chemical potential of 
the right lead (Fig. 6 d) constitutes a measurement of 
that superposition (collapse of the wave- function) : the 
electron is either in the right dot (with a high probabil- 
ity 1 — ^S{az)f) and tunnels out, or the electron is in 
the left dot (and nothing happens because the left level 
is still below fj, and the system is Coulomb blocked). 

If the tunnel rates Tji^Tl to the right and left leads 
are sufficiently larger than the inverse of the cyle duration 

^cycle, 

rfl,ri»i-ii^ (33) 

the decharging of the right dot and the re-charging of 
the left dot from the left lead is fast enough to bring 
the system back into its initial state with one additional 
electron on the left dot. In this case, the precise value of 
TfijTL, and the precise shape of the £(t)-pulse for < 
t < ^cycie has no effect on the total charge transfered 
within one cycle. Then, since the probability to transfer 
one electron from the left to the right in one cycle is given 



by 1 — ^S{az) f, on the average an electron current 

(/> = -ei^i^ (34) 

'^cyclc 

flows from left to right. Note that the leads essentially 
act as classical measurement devices of the quantum- 
mechanical time-evolution between the two dots. Mea- 
suring the current (/) as a function of the pulse length 
tf — n/n then offers a scheme to make quantum mechan- 
ical oscillations such as those predicted in Eq. (31) visible 
in the electronic current, similar to the recent experiment 
by Nakamura et al. in a superconducting Cooper pair 
box^s. 



D. Discussion 

The 'Rabi- pulse', Eq. (27), keeps the energy difference 
to the excited state |-|-) constant throughout the adia- 
batic rotation. In this case, it follows that the dissipative 
contribution to 6{az)f, Eq.(31), is due to p{u) = A/fi), 
i.e., one fixed boson frequency only. Furthermore, the an- 
alytic result, Eq.(31), in principle allows to extract the 
value of the phonon spectral density, p(A), for 'Rabi'- 
pulses at fixed energy difference A. 

Another important observation is the fact that the ro- 
tation remains dissipation-free if A is chosen to coincide 
with a zero of p{uj). This defines a 'decoherence-free man- 
ifold' in the parameter space of the system. Zeroes in 
p{uj) at certain frequencies u — ujq have been predicted'^^ 
in free-standing, two-dimensional phonon cavities (slabs). 
In fact, the control of vibrational properties of quan- 
tum dot qubits has been suggested*^, and considerable 
progress has been made in the fabrication of nano- 
structures that are only partly suspended or even free- 
standing^^"^^. 
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In a thin- plate cavity modcP^, the zeroes in p{uj) are 
due to symmetry and geometrical confinement both for 
deformation potential and piezoacoustic phonon scatter- 
ing in second order of the coupling constant g. This 
means that within a phonon cavity, the dissipative con- 
tribution to 6{az)f, Eq. (31), can be 'switched off' if A 
is tuned to the energy fuvo of a decoupled phonon mode. 

For GaAs slabs of width 1/xm, this energy has been 
predicted to be of the order Hluq w lOfieY. An energy 
A = 10/LteV corresponds to ksT « lOOmK such that 
the temperatures shown in Fig. 4 are attainable in cur- 
rent experiments. A typical pulse frequency Q/2'k = 
0.01A/27r/i then would correspond to 24 MHz, the corre- 
sponding current being of the order (/) ^ eil/27r « 4pA. 
Furthermore, we mention that the weak coupling regime 
with g « 0.01 ~ 0.05 seems to be justified to describe re- 
cent experiments on phonon coupling in double quantum 
dots3>5. 



IV. SPIN-QUBIT SWAP AND DECOHERENCE 
IN TWO-QUBITS 

In this section, we apply our formalism to spin-qubit 
swaps and charge decoherence in two-qubit operations. 
These requires larger Hilbert spaces and can be realized 
by two electrons on a double dot. Loss and DiVincenzo 
have introduced a detailed scheme for 'quantum compu- 
tation' with spin-states of coupled single-electron quan- 
tum dots®. Dephasing of spin degrees of freedom due 
to spin-orbit coupling or the coupling to nuclear spins is 
expected to be much weaker than dephasing of charge de- 
grees of freedom. Nevertheless, spin and charge become 
coupled during switching operations whereby charge de- 
phasing also effects spin-based qubits^'''. 

Here, we consider a specific two-qubit swap operation 
as discussed recently by Schliemann, Loss and MacDon- 
ald: two electrons with spin are localized on two coupled 
quantum dots A and B, giving rise to a basis of six states. 
During the operation charge decoherence occurs for in- 
termediate states that get involved in the swap operation 
when charge is tunneling between the dots. Piezoelectric 
phonons then couple to the electron charge and incoher- 
ently mix states in the singlet sector which leads to a 
loss of fidelity of the swap operation. As in the previous 
section, we only consider the coupling of phonons to the 
charge density which yields a lower bound for dephasing. 

A. Two-Qubit Hamiltonian 

The four basis vectors with the two electrons on differ- 
ent dots are the spin singlet and triplets 

l^i) := 2-i/2(ct^^4^_c^^4^)|0) 
|r-i) := cl^cl^lO), |Ti) :=c^^4^|0) 
|T°) := 2-V2(,t ,t +J ,t )|o). (35) 



The remaining two states with two electrons on dot A 

('left') or dot B ('right') are 

\L) ■■= c\^c\^\0)=2-'/mS2) + \S,)] 

\R) := 4^c;j0) = 2-1/2 [|52)-|53)] (3g) 

which are superpositions of two spin singlets |<S'2,3), 

|52,3) := 2-V2(ct^^ci^± 4^4^)10) (37) 

that differ in their orbital wave function. 

As in^^, wc specify the two-qubit swap as an adiabatic 
rotation from an initial state \i) to a final state |/), 

\i) := i= [|T0) + \S,)] - I/) := i= [IT") - \S,)] (38) 

that can be achieved®'^ by an adiabatically opening and 
then closing of the tunnel barrier between the two dots as 

a function of time. This operation leads out of the sub- 
space span{|S'i), |T°)} since it involves intermediate dou- 
bly occupied states in span{|L), \ R)} (= span{|52), jS's)}) 
which altogether define a four-dimensional Hilbert space, 
Ti.^^\ in which the two-qubit swap takes place. To be 
more specific , we use a basis of n'-^^ as given by the 
three singlets \Sj) and the triplet |T°), 

|0):=|T0), \j):=\Sj), (i = 1,2,3), (39) 

to define the time-dependent Hamiltonian 

= j2'j\j){j\+m)[\i){2\ + \2){i\], m 

i=o 

where Sj denotes the energies of the spin singlet states, 
ei = eo, S2 = so + Uh, £3 = eo + Uh — 2X with the spin 
triplet energy eo, the on-site Coulomb repulsion Uh > 0, 
the exchange term X > 0, and the time-dependent tunnel 
coupling element between the dots Tc{t). 

B. Electron-Boson Coupling, Master Equation 

The total Hamiltonian in presence of bosons coupling 
to the charge degree of freedom is derived in Appendix 

B, 

H^^\t) = H^^\t) + ^<7,i + Hb, (41) 

it has exactly the same form as in the one-qubit case, 
Eq.(6), but with the free Hamiltonian HQ^\t) replaced by 

(2) 

-fffj (t), the coupling constants gQ replaced by gq, and 
az \L){L\ ~ \R){R\ now referring to the two-particle 
states Eq. (36). In analogy with the one-qubit case, 
we introduce the matrix elements of the reduced density 
operator in the interaction picture with respect to Hp by 

Pi,(i) :=lVph(j|p(t)|i), i = 0,l,L,i?. (42) 
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The coherent part of the time evolution of the Pij{t) is 
obtained trivially from the Liouville equation, using the 

Hamiltonian H^\t), Eq. (40). The non-trivial part are 

the additional terms due to the electron-phonon coupling. 
A systematic perturbation theory in the latter starts from 
the four adiabatic eigenstates of the unperturbed Hamil- 
tonian H^^\t), Eq. (40). 

Here, we restrict ourselves to the case of small tunnel 
coupling, 



\Tc{t)\ < Uh.2X, 



(43) 



where inelastic transitions are determined by the dynam- 
ics in the subspace spanned by the states |2) and |3) and 
admixtures from |1) through the hybridization between 
|1) and |2) can be neglected. This case is particular- 
ily interesting since the (adiabatic) energy difference 2X 
between |2) and |3) remains constant throughout the op- 
eration. As we show below, if 2X coincides with a zero 
hwi) of the boson spectral density, the operation is again 
dissipation- less in second order of the electron-boson cou- 
pling. 

Within the Born and Markov approximation one finds 




200 300 
timet 

FIG. 7: Expectation value 1 — {f\p{t)\f) for a two-electron 
double dot qubit swap from initial state \i) to final state ]/), 
Eq. (38), as a function of time (in units of Un/h, Uh' on- 
site Coulomb repulsion) . The corresponding dcphasing times 
{g > 0) are (from top to bottom) F'^ = 84, F'^ = 169 
and = 635. Inset: (i|p(t)|i) and {f\p{t)\f) for vanishing 
phonon coupling g = 0. 



PLfllep = -TpLR-T-PRR + T+pLL 

PLilep = -^PLz- —Pm, i = 0,l 

_ r 

Pmlep — 



0,1. (44) 



Here, the rates F and T± are defined in analogy with 
Eq. (14), 



r = 27r^p(A)coth(^^ 

r. = 4fp(A), 



(45) 



with the energy difference A given by A = 2X. Notice 
that due to the limit Eq. (43), A does not depend on 
time. 

In Fig. (7), we show results of a numerical evaluation 
of the two-electron master equation. The spectral density 
Ppz{u}), Eq. (7), was used to obtain the rates, Eq. (45), for 
different coupling strengths g and energies hiOd {oJd = c/d, 
with c the speed of sound and d the dot center distance). 
Using Uh = 1 as energy scale, our results coincide with 
those of Schliemann et al. for 5 = with a pulse of the 
form^^ 



Tc{t) = 



1 + cosh(t/T)/ cosh(T/2T) 



(46) 



as shown in the inset. Here, we chose To = 0.05, T = 400, 

T = 50, X = 0.5, and the temperature 1/(3 — 0.1. A 
value of Uh = ImeV corresponds to a time of 4 x 10~^^ 
s. 

The diagonal element {,f\p{t)\,f) is a measure of the 
fidelity of the swap operation, transfering \i) into |/), 



cf. Eq. (38). Even in absence of dissipation, the non- 
adiabacity of the operation results in a finite value of 
1 — {f\p{t)\f) after the swap^^. The electron-phonon in- 
teraction acts when charge between the dots is moved 
during the opening of the tunneling barrier. Conse- 
quently, the two states |2) and |3) become mixed inco- 
herently, leading to a finite, irreversible occupation prob- 
ablity of the energetic lower state |3) even after the pulse 
operation. We mention that the data shown here corre- 
spond to a worst case scenario where the inelastic rate 
r, Eq. (45), is essentially determined by a large en- 
ergy difference 2X between |2) and |3) which enters into 
p(A = 2X). Then, spontaneous emission of phonons oc- 
curing during the slow swap (cf. Eq. (43)) leads to a 
dephasing rate F « ngX. In this case, even relatively 
small values of 3 (5 < 0.05 in experiments on lateral 
dots-^'^^) can lead to a considerable fidelity loss of the 
operation. 

On the other hand, our results also demonstrate that a 

complete suppression (at least to lowest order in g), simi- 
lar to the one-qubit case discussed above, is theoretically 
possible for the two-qubit swap in phonon cavities. If A 
is tuned to the energy Hloq of a decoupled phonon mode, 
the inelastic rates, Eq. (45), can be 'switched off and 
the operation again is 'dissipationless'. 



V. CONCLUSION 

In conclusion, we have studied adiabatic quantum state 
rotations in coupled quantum dots in presence of dissi- 
pation. For the one-qubit case, we suggested a scheme 
to extract quantum oscillations in the inversion change 
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5{iJz) f from the average pump current through the dots, 
Eq. (34). The analytic result, Eq.(31), furthermore al- 
lows to extract the value of the phonon spectral density, 
/9(A), for 'Rabi'-pulses at fixed energy difference A. 

For the two-qubit case, we have quantified the effect 
of dissipation on a swap operation, where states in the 
singlet sector are incoherently mixed. Although we have 
only considered the case of small tunnel coupling between 
the dots, as in the onc-qubit case wo have found that two- 
qubit rotations can be performed such that dissipation is 
due to p{u) = A/h), i.e., one fixed boson frequency only. 
If A is chosen to coincide with a zero of p{lli) as in free- 
standing phonon cavities, this defines a 'decoherence-free 
manifold' in the parameter space of the system. 
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0144, the German DFG project Br 1528/4-1, the UK 
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cuits Network. Discussions with R. H. Blick, S. De 
Franccschi, E. M. Hohberger, L. P. Kouwenhoven, and 
F. Renzoni are gratefully acknowledged. 



APPENDIX A: CHARGE-PHONON COUPLING: 
ONE-QUBIT 

Phonons coupling to the electron charge in a double 
quantum dot are described by a spin-conserving electron- 
phonon Hamiltonian 



Hp, 



qPq 



a_Q -I- a, 



(Al) 



where Uq is the creation operator for a phonon mode 
Q, po the Fourier transform of the operator /5(x) = 
of the electron density, and Aq = X'^q 
is the microscopic electron-phonon interaction matrix el- 
ement. 

The electron-phonon interaction potential in real 
space in first quantization is given by Kip(x) = 

^qs^^^ ("^-Q + "^q) • consider the simplest case 
of an electron-phonon interaction with only diagonal 
terms 



Hep = ["Q^^ + "Q^« 



Q 



where iV, = \i){i\ {i = L, R), a^^ = J^a- "^Qo"' ^^'^ 



(A2) 



L/R 



\q J d^^e^'^-pL/nA^). (A3) 



Here, pi,<j(x) := {i\'^l{x)'i'crix) i ~ L,R, is the den- 
sity of electrons with spin a in the left (right) dot as 
obtained from the (many-body) ground states \L) {\R)). 

Assuming two identical dots with pi<j(x) = 
Pct(x — r^/zj) to be smooth functions centered around 



the left (right) dot centers r^/ji, we obtain 

«q'' = ^Qexp(iQri/;j)Pe(Q), (A4) 

where Pe(Q) = Y^a I d^^e^^^Pai'^) is the 'form factor' 
of the electron density in the left and in the right dot. 
In the electron-phonon interaction, Eq. (A2), we write 

"Q^i + "Q^fl = I Wzia^ - a§) + i(a^ + a§)](A5) 

and recognize that the term proportional to the unity 
operator 1 in the two-dimensional Hilbert space H^^^ 

leads to a mere renormalization of the total energy by 
the phonons in Eq. (A2): the energies in both dots 
are shifted by the same amount. In fact, if the coupling 
constants and ctq were the same, Eq.(A2) would de- 
scribe a trivial coupling of the phonons to Nl + Nn = 1 
which would not affect the dynamics of the system, while 
here (due to Eq. (A4)) we find the phase relation 



(A6) 



where d is the vector pointing from the center of the left 
to the center of the right dot. 

Defining qq = ctQ — aQ, the term ^<7z{ctQ — ctQ) in Eq. 
(A5) yields the interaction term ^cr^ A in the Hamiltonian 
J?(i)(i), Eq. (6). 



APPENDIX B: THE RATES 7 AND -y±, EQ. (14) 

The ohmic form Eq.(8) with s = 1 for the spectral 
function p[uj) = gu)e~'^/'^'= is used in the discussion of 
the one-qubit dynamics in section II. In this case, the 
corresponding rates 7 and 7±, Eq. (14), are given by 

. = fl(f! + 2r^A.-M»..(f) 



Re{7±} = g'^ ( 4 - I A e-^'^^ coth 



A2 V/3 2 



2 



Im{7±} = =F5 



eTcU>c 



Tc r 

2A2 7o 



du> (je "^"^ 



{ e/?'^ - 1 ( 



2A2 

At£ A±£ 



e'^'^ A oj + A 

1 /A±£ A^e 



\ui-A oj + A 



(Bl) 



Note that the ohmic case s = 1 is peculiar in the sense 

that a contribution cx 1//3 = T linear in the temper- 
ature T appears in 7 and Re{7±} due to the limit 
lim^^ p{u})nB{u)) under the integral in Eq. (14), which 
gives a finite contribution for s = 1. However, the mi- 
croscopic calculation^" of ppz(w), Eq.(7), shows that the 
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ohmic spectral density as an 'envelope' of pp^iuj) fails to 
be correct for small uj where Ppz{io ^ 0) cx cj"^ for piezo- 
electric phonons in double quantum dots. The contribu- 
tion to 7 and Re{7±} linear in T therefore is unphysical 
but of minor relevance in the numerical calculations any- 
way, as we have checked. 



2"^/^ [152) ± IS'a)] with two electrons on the left (right) 
dot, Eq. (36). Using the completeness relation in Ti^^^ 
yields 

("Q + «q)[|0)(0| + |1)(1| + 2a^\L){L\ + 2a^\R){R\ 
= i(a^+a^) + (a^-a^)a„ (C4) 



APPENDIX C: CHARGE-PHONON COUPLING: 
TWO-QUBIT 



where again we defined cr^ \L){L\ — \R){R\. 
The total Hamiltonian now can be written as 



As in the one-qubit case, we start from an electron- 
phonon interaction Hamiltonian with only diagonal terms 



= J2 "Q-S>Sa 

QrT;j=A,B 



a_Q + a, 



(CI) 



where the matrix element 

< ■■= Aq / rf3xe^Q-01*t(x)*,(x) \j) (C2) 

is calculated with the single particle states \j = A,B). 
We assume that the single electron densities in Eq.(C2) 
are spin-independent, ccq^ = Uq whence 



Q 



o_Q -I- 



(C3) 



with the number operators Nj := ^ja^jcr- '^^^ rion- 
vanishing matrix elements of the number operators are 
{j\Na\j) = l,mB\j) = 1 where j = 0, 1, {L\Na\L) = 2, 
and {R\Nb\R) = 2. Here, we used the states \L/R) := 



H^''\t) = Y. e^\3){j\ + ^[\l){{L\ + {R\) + H.c] 

3=0,1, L,R ^ 



(C5) 



where A = J2q 9q. (o-q + «q) , El/k = ^ (^2 + £3) and 
(Tx = \L){R\ + |-R)(-^|. Comparing the coupling constant 
9Q = ctq— ttq in the one-qubit case and = 2(aQ— ag) 
in the two-qubit case, and keeping in mind the definitions 

for Q.Q^ and oc^^ , Eq. (A3) and Eq. (C2), we realize 
that for spin- independent electron densities fulfilling 

(L|vl/t(x)vl/,(x) \L) = (A|vl/t(x)vl/,(x) 1^) 

(i?|vE't(x)*,(x) \R) = (S|vl/t(x)vl/,(x) \B) (C6) 
we obtain identical coupling constants gQ = 5q: the spin 
sum in Eq. (A3) yields a factor 2. An exception is the 
one-qubit with only one electron having a fixed spin (so 
that there is no spin sum). In that case and with Eq. 
(C6) holding, the one-qubit coupling is half the two- 
qubit coupling qq. 
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